Symmetry breaking in dipolar matter-wave solitons in dual-core couplers 
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We study effects of the spontaneous symmetry-breaking (SSB) in solitons built of the dipolar Bose- 
Einstein condensate (BEC), trapped in a dual-core system with the dipole-dipole interactions (DDIs) 
and hopping between the cores. Two realizations of such a matter- wave coupler are introduced, 
weakly- and strongly-coupled. The former one in based on two parallel pipe-shaped traps, while the 
latter one is represented by a single pipe sliced by an external field into parallel layers. The dipoles 
are oriented along axes of the pipes. In these systems, the dual-core solitons feature the SSB of the 
supercritical type and subcritical types, respectively. Stability regions are identified for symmetric 
and asymmetric solitons, and, in addition, for non-bifurcating antisymmetric ones, as well as for 
symmetric flat states, which may also be stable in the strongly-coupled system, due to competition 
between the attractive and repulsive intra- and inter-core DDIs. Effects of the contact interactions 
are considered too. Collisions between moving asymmetric solitons in the weakly- symmetric system 
feature elastic rebound, merger into a single breather, and passage accompanied by excitation of 
intrinsic vibrations of the solitons, for small, intermediate, and large collision velocities, respectively. 
A T^T-symmetric version of the weakly- coupled system is briefly considered too, which may be 
relevant for matter-wave lasers. Stability boundaries for T^T-symmetric and antisymmetric solitons 
are identified. 

PACS numbers: 42.65.Tg; 03.75.Lm; 47.20. Ky; 05.45.Yv 



I. INTRODUCTION 

Studies of Bose-Einstein condensates (BECs) made of dipolar atoms or molecules had produced a great deal of 
fascinating experimental and theoretical results, which were summarized in recent reviews [H,!!! and [2|,[3[, respectively. 
The continuation of the work in this direction has yielded new remarkable findings, such as the prediction of various 
pattern- formation mechanisms 0, HI (which share some features with the formation of patterns in ferrofluids [6[), 
analysis of the stability of the dipolar BEC trapped in optical-lattice (OL) potentials [?! and of the roton instability 
[s,, the possibility of the Einstein - de Haas effect [9], etc. Important experimental achievements, which offer 
new perspectives for studies of dipole-dipole interactions (DDIs) in atomic condensates, are the creation of BEC in 
dysprosium [Io| and erbium [11]. Parallel to that, essential results have been obtained for degenerate quantum gases 
of dipolar fermions [s, 12]. 

In addition to their own physical significance, dipolar condensates may also be used as quantum simulators [l3|, [l^ 
representing other physical media where nonlocal nonlinearities play a fundamental role. These include the heating 
and ionization of plasmas [l5| , nonlinear optics of nematic liquid crystals [l6| , of waveguides sensitive to temperature 
variations [l^, and of semiconductor cavities [18], BECs with long-range interactions induced by laser illumination 
[l9|, and others settings. 

An interesting ramification of the study of collective nonlinear modes in the dipolar BEC is the prediction of solitons 
(which have not yet been reported in experimental works). In effectively one-dimensional (ID) traps, solitons were 
analyzed in both continual [20] and discrete ^] settings, the latter one corresponding to the fragmentation of the 
BEC by a deep OL. In a similar form, ID solitons supported by the attractive DDIs were predicted in the Tonks- 
Girardeau gas of dipolar hard-core bosons [22|. Taking into account the 3D structure of the quasi-lD cigar-shaped 
traps, the solitons, including ones with embedded vorticity (cf. similar modes introduced earlier in the context of 
BEC with local interactions |23|) were further studied in Ref. [IJ, and gap solitons in a similar setting, but including 
an OL potential, were considered too psj. In the 2D system, discrete fundamental solitons and solitary vortices with 
long-range DDIs between sites of the lattice can be constructed easily [26|. In the continual 2D model, fundamental 
p7| and vortical [28] solitons were constructed in the isotropic setting, assuming that the sign of the DDI could be 
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reversed from repulsion to attraction by means of a rapidly oscillating ac field [29|. 2D fundamental and vortex 
solitons supported by a trapping potential were introduced in Ref. [30]. 

Without reversing the DDI sign, stable 2D anisotropic solitons, corresponding to the in-plane polarization of dipoles, 
were constructed in Ref. [31], by means of the variational approximation and systematical numerical simulations. The 
variational approximation for 2D solitons supported by the DDI was analyzed in Ref. ^2| , and a rigorous proof of the 
existence of such solitons was provided too Also studied were more complex situations, such as the formation of 
a multi-soliton patterns as a result of the development of the modulational instability of an extended state [33^. 

The long-range character of the DDI makes it possible to consider interactions between condensate layers trapped 
in parallel planar waveguides. The DDI couples them by nonlinear forces even in the absence of hopping (tunneling) 
of atoms across gaps separating the layers [5| (the isolation of parallel layers can be provided by a strong OL field 
whose axis is perpendicular to the layers [35|)- This nonlocal interaction gives rise to "indirect" scattering of 2D 
solitons moving in the separated layers [36[ , and the formation of bound states of such solitons [s^ • The creation of 
multi-soliton filaments and checkerboard crystals in multi-layered stack was predicted too [s^. 

The model of nonlocal DDIs between parallel layers considered in Refs. [HI, IsgMsS^ did not take into regard the 
hopping (tunneling of atoms, alias linear coupling) between the layers. On the other hand, models of dual-core 
couplers^ with intrinsic local nonlinearity acting in both cores, were studied in detail in terms of optics and matter 
waves, starting from the analysis of the spontaneous symmetry breaking (SSB) of CW (continuous- wave, i.e., uniform) 
states in dual-core optical fibers with the cubic and more general forms of the intra-core nonlinearity. In that system, 
the linear coupling is caused by the overlap of the evanescent field, originating from each core, with the parallel one. 
The SSB happens in the dual-core fiber, as a result of the interplay of the linear coupling and intrinsic nonlinearity, 
with the increase of the total power of the CW beam. The analysis of the SSB was extended, in full detail, to 
temporal and spatial solitons |4ll , |42| . and to optical domain walls [43|- The SSB effects were also studied for solitons 
in dual-core fiber Bragg gratings [44^, in two-tier waveguidingarrays (for discrete solitons) [45[, in parallel-coupled 
waveguides with the quadratic (second- harmonic-generating) |4^ and cubic-quintic (CO) [l^ nonlinearities, as well as 
for dissipative solitons in linearly-coupled CQ complex Ginzburg-Landau equations [43 • Recently, a similar analysis 
was developed for the SSB of solitons in T^T-symmetric couplers, with mutually balanced loss and gain (and identical 
cubic nonlinearities) acting in the two cores |50i]. Unlike the above-mentioned settings, in the latter case the SSB 
destroys symmetric solitons, rather than replacing them by stable asymmetric ones. The linear coupling in optics may 
also represent the mutual interconversion of two polarizations of light in twisted fibers (which, in particular, are used in 
the so-called rocking filters) [5l|, twisted photonic-crystal fibers [5^], twisted fiber gratings [53|, or the interconversion 
of two waves with different carrier frequencies, caused by the electromagnetically-induced transparency [5^ . 

Similar dual-core (ahas double-well) settings, approximated by linearly-coupled Gross-Pitaevskii equations (GPEs), 
were introduced for the mean- field wave functions describing trapped BEG with local interactions [55| . A similar linear 
coupling accounts for the mutual interconversion in a mixture of two different atomic states, induced by a resonant 
electromagnetic wave [56] . The linearly coupled GPEs were used to predict the shift of the miscibility-immiscibility 
transition in BEG or fermionic mixtures of two states connected by the linear interconversion [57] and the stabilization 
of 2D solitons against the collapse in a linearly-coupled binary system with attractive and repulsive intra-component 
interactions [58]. The SSB of matter- wave solitons trapped in ID and 2D linearly-coupled cores was analyzed in Ref. 
[55[. In the ID situation, a more accurate description, which, nevertheless, yields similar results for the solitons' SSB, 
is provided by the two-dimensional GPE, which, instead of postulating two ID wave functions in the two parallel 
cores with the linear exchange between them, introduces a single 2D wave function comprising both cores [5^]. This 
includes the case when the cores are defined by nonlinear pseudopotentials (rather than by a linear trapping potential), 
i.e., local modulation of the self-attraction coefficient |6Qj . 

Although the previous works analyzed many aspects of the SSB in dual-core GW and solitonic states, those works 
were dealing solely with local intrinsic nonlinearities. Only in a very recent paper [6l|, a shift of the SSB transition of 
solitons in the coupler with nonlocal nonlinearity of the thermal type, typical to optical systems |17|], was considered. 
The analysis was performed for two opposite Hmit cases, viz., the weak nonlocality characterized by a small correlation 
radius, which may be approximated by the first two terms of the expansion of the nonlocal cubic term, and the opposite 
limit of the infinite correlation radius, which corresponds to the two-component quasi-linear model of "accessible 
solitons" [H. 

The aim of the present work is to consider the SSB of solitons in the effectively ID dual-core coupler filled by the 
dipolar condensate, which exhibits the interplay of the long-range DDIs and linear hopping between the quasi- ID 
cores. In particular, the DDIs act both inside the cores and between them, while the thermal nonlocality considered 
in Ref. [6i| could not act across the gap separating the parallel waveguides. It is relevant to stress that, in the 
absence of the longitudinal dimension, the double- well setting is not sufficient to exhibit the nonlocal character of the 
interactions in the dipolar condensate, the minimum necessary configuration being based on a set of three potential 
wells [HI. 

Two coupler configurations are considered here, as shown in Fig. [TJ The first setting, presented in Fig. l(a)[ is 
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FIG. 1: (Color online) (a) In the weakly-coupled system, the dipolar condensates are trapped in parallel pipes of diameter 6, 
separated by distance a. The arrows represent the orientation of the dipoles. (b) The strongly- coupled system is shown by 
means of the cross-section image of the condensate trapped in the single pipe of diameter 6, which is sliced by a repelling laser 
sheet into two layers, with effective separation a between them. Symbol represents the orientation of the dipoles, which are 
perpendicular to the figure's plane. 



based on two identical condensate-trapping pipes of diameter 6, separated by distance a. In this case, b < a i s implied, 
hence the system may be naturally called a weakly- coupled one. The other setting is shown in Fig. |2(b)[ with the 
condensate loaded into the single pipe of diameter 6, which is sliced into two parallel layers by a thin potential barrier 
of small thickness a. The barrier can be induced by a repulsive light sheet (blue-shifted one, with respect to the 
atoms) [64]. The latter setting implies 6 > a, and it will be named, accordingly, a strongly- coupled system. In either 
case, the dipoles are polarized along the pipes' axes, hence the DDIs are attractive inside the cores, which makes it 
possible to form solitons in each one f20|. 

Aiming to study the SSB, alias the symmetry-breaking phase transition, in these settings^ is relevant to recall that 
there are two types of the SSB bifurcation, namely, the subcritical and supercritical ones [65|, which are tantamount 
to the phase transitions of the first and second kinds, respectively. In the subcritical situation, branches representing 
asymmetric modes emerge as unstable states at the bifurcation point, then go backward in the bifurcation diagram, 
and get stabilized after turning forward. In the supercritical setting, the asymmetric branches emerge as stable ones 
at the bifurcation point and immediately continue forward. 

The rest of the paper is organized as follows. In Sec. II, we construct symmetric solitons in the two settings 
presented in Fig. 1, and then identify the SSB transitions to asymmetric solitons. Antisymmetric solitons and their 
stability are considered too (these solitons do not undergo any bifurcation). Basic results for the stability of different 
modes are presented (including a flat state, which may also be stable in the strongly-coupled system). Effects of the 
local (contact) nonlinearity on the SSB are considered too. In Sec. Ill, we study collisions between moving asymmetric 
solitons. In Sec. IV, a T^T-symmetric [66| extension of this system is briefly considered, which includes gain and 
loss applied to the two cores (similar to the T^T-symmetric coupler with the local nonlinearity introduced in Refs. 
HH^IEqI)- The paper is concluded by section V. 



II. THE SYMMETRY-BREAKING BIFURCATION AND STABILITY OF SOLITONS 



A. The coupled Gross-Pitaevskii equations 

In the usual mean- field approximation [H, Q, both dual-core settings introduced in Fig. [1] are described by the 
system of ID linearly-coupled GPEs for the wave functions in the two cores, ipi and ?/^2- In the scaled form, the 
equations are 
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x'P)^/^, see Fig. [TJa). In this notation, is the coupling parameter 



(hopping coefficient), g represent the local interaction (repulsive in the case of ^ > 0), and the orientation of the 
dipoles in Fig. [T] corresponds to Gdd > 0, which implies the attraction between dipoles in the given core and 
repulsion between the cores for cos^ < 1/3. The first term in the integrand, that accounts for the intra-core DDIs, 
is regularized by the transverse diameter, 6, which is an approximation sufficient for producing ID solitons (2Qj . 
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The solitons will be characterized by the total norm (proportional to the number of atoms in the condensate), 
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To focus on SSB effects dominated by the DDI, we will first drop the local nonlinearity, setting ^ = (in the 
experiment, this can be done by means of the Feshbach resonance [S^]); effects of the contact interactions will be 
considered afterwards. Then, we scale the units to fix = 1 and Gdd = 1, the remaining free parameters being a, 6, 
and P (and g too, in the end). 

Stationary solutions to Eq. ([1]) with chemical potential /i are looked for in the usual form, '^1,2 (^7^) = 
exp (— i/it) 0i^2(^), with real functions 0i,2(^)- In particular, ^i,2(^) = for symmetric solutions obeys the 

stationary equation. 
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where Hi = Gdd = 1 is fixed, as said above. Equations ([T]) and (|2j) are solved below by means of numerical methods. 
In particular, stationary solutions were found below by means of the imaginary-time propagation method [68| with 
periodic boundary condition, while the stability of the solutions was tested by means of the integration in real time. 



B. The weakly-coupled system (a > b) 

Typical examples of stable symmetric, asymmetric, and antisymmetric solitons found in the system represented 



by Fig. |l(a) are displayed in Fig. [2j Naturally, the asymmetric solitons, morphed by the stronger nonlinearity, are 
narrower and taller than their symmetric counterparts. 

It is relevant to mention that, unlike the standard model of the coupler with the local cubic nonlinearity, where 
symmetric and antisymmetric solitons are available in an obvious exact form, and asymmetric ones can be effectively 
described by means of the variational approximation [69| , in the present strongly nonlocal system such an approx- 
imation is not practically possible. Nevertheless, some results can be obtained in an approximate analytical form for 
the present system too, see Eqs. ([5]) and (|6]) below. 

Before proceeding to the consideration of the transition between the symmetric and asymmetric solitons, we display 
the stability area of the antisymmetric ones (which do not undergo any antisymmetry-breaking bifurcation), in the 



plane of (h^P)^ in Fig. 2(c)[ The stability boundary may be fitted to curve 

P = 7b^. (3) 

Unstable antisymmetric solitons gradually decay into radiation (not shown here in detail). 

The SSB for solitons is summarized in Fig. [2] by means of the bifurcation diagrams, which clearly show that 
the symmetry-breaking bifurcation, driven by the nonlocal attractive DDIs, is supercritical, on the contrary to the 
commonly known subcritical bifurcation in the coupler with the local self- focusing nonlinearity [4l|, |42[ . A similar 
trend to the change of the character of the SSB bifurcation for solitons from sub- to supercritical, with the increase of 
the degree of the nonlocality, was recently demonstrated in the model of the optical coupler with the weakly nonlocal 
thermal nonlinearity [6l|. 

As seen from the structure of the first term in the integrand of Eq. ([T]), the decrease of diameter b of the parallel 
pipes implies effective enhancement of the nonlinearity. This, as well as the direct strengthening of the nonlinearity 
due to the increase of the total norm (P), leads to the symmetry breaking, as seen in Figs. [3(b)] and [3(a)| Further, 



Fig. 3(c) I demonstrates the related effect of the decrease of the critical value. Per, of the total norm at the bifurcation 



point with the decrease of b. The latter dependence may be fitted to formula Per = 36^, which is explained by the fact 
that, at small 6, the nearly diverging first integral term in Eq. ([2]) may be estimated as where A is the soliton's 

amplitude. Its balance with other terms in the equation leads to estimates for scalings of the amplitude and width: 

P/b,W ^b^/P. (4) 

Then, as the SSB point is determined by the competition between the nonlinear intra-core terms and linear inter- 
core-coupling ones in Eq. ([T]) [4^-[4l|, the corresponding scaling for the value of P at the critical point indeed takes 
the form of 

Per - V^b^ (5) 
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FIG. 2: (Color online) Examples of stable symmetric (a), asymmetric (b), and antisymmetric (c) solitons found in the weekly- 
coupled system [see Fig. [B^a)] for a = 1, 6 = 0.4, and total norm P = 0.4 (a) and P = 0.6 (b,c) The antisymmetric solitons are 
stable below the boundary, P ^ 76^, shown in the plane of (6, P) in panel (d), at different values of a. 



(here, k is kept for clarity, although it was actually scaled to be z^: = 1). 

Note that, according to Eq. (|3]), the stability region for antisymmetric solitons is much smaller at small b. This 
agrees with the general trend of the antisymmetric solitons in nonlinear couplers to be more fragile modes than their 
symmetric counterparts [4l|, Hil, [50| , due to the obvious fact that they correspond to a larger coupling energy. In 
fact, the cubic scaling in Eq. (|3]) may be qualitatively explained too, although in a more vague form than Eq. (|5]). 
Indeed, the antisymmetric soliton with amplitude A is subject to an oscillatory instability characterized by complex 
growth rates (eigenvalues), A = ±m + Re (A), where, in the generic case, an estimate Re (A) ~ A is valid (see, e.g. 
Ref. [50|). Because the instability is oscillatory, the corresponding perturbations tend to escape from the region of 
width W , occupied by the soliton [see Eq. dU], within time r ~ W/Vgr^ where the group velocity is determined by 
the characteristic wavenumber of the perturbation mode, k ~ W~'^^ i.e., r ~ W^. The instability accounted for by 
the escaping perturbations seems as convective one, that will have enough time to destroy the antisymmetric soliton 
under condition Re (A) r ~ 1, i.e., according to the above estimates, 

AW^ ^b^/P ^1, (6) 

which qualitatively explains the numerically found fit ([3j). 

It is also worthy to note that additional analysis demonstrates that, as is strongly suggested by Figs. |2(d)| and 3(c)[ 



the general picture of the SSB, being sensitive to the value of 6, shows little dependence on the separation between 
the cores, a (in the case of a > 6, which is considered here). In other words, the SSB in the weakly-coupled system, 
in accordance with its name, is weakly sensitive to the DDI between the parallel cores, which renders the picture 
relatively simple. 
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FIG. 3: (Color online) (a) The bifurcation diagram for solitons in the weakly- coupled system: The soliton's asymmetry, measured 
by the deviation of the share of the total power in one core (Pi/P) from 0.5, versus total norm P. (b) The bifurcation diagram 
as a function of the pipes' diameter, b. The circles located along the solid and dashed lines represent stable and unstable 
solutions, respectively, (c) The critical value of the total norm at the symmetry-breaking point as a function of b, for fixed 
a = 1. Stable asymmetric solitons exist above the curve, while the symmetric solitons are stable below it. The curve is well 
fitted by Per = 36^ 
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FIG. 4: (Color online) Examples of stable modes found in the strongly-coupled system, (a) A symmetric soliton , with 
parameters (a, 6, P) — (0.3, 0.5, 0.2). (b) A symmetric flat state with (a, 6, P) — (0.1, 0.5, 0.2). (c) An asymmetric soliton, with 
(a,6,P) = (0.2,0.5,1). 



C. The strongly-coupled system (a < b) 



In the setting displayed in Fig. l(b)[ the small separation between the effective cores makes effects of the inter- 



core DDIs essentially stronger, in comparison with the weakly-coupled system. In fact, the strongly-coupled system 
realizes an example of competing interactions^ namely, intra-core attraction and inter-core repulsion. A somewhat 
similar example is the discrete Salerno model with competing signs of the onsite and intersite cubic nonlinear it ies, 
which was studied in ID and 2D settings [70]. 

The numerical solution demonstrates that, in addition to symmetric and asymmetric solitons, stable flat states 
with unbroken symmetry between the cores also exist in the strongly-coupled system, being stabilized by the strong 
dipolar repulsion between the cores. A typical example of such stable states is displayed in Fig. |4l 

In comparison to their counterparts in the weakly coupled system (cf. Fig. [2]), the symmetric solitons are wider in 
the case of the strong coupling. This is also an effect of the repulsive DDI between the cores, which partly cancels the 
intra-core attraction, which forms the symmetric modes. The transition to the asymmetric soliton occurs, as before, 
with the increase of the total norm, the difference from the weakly coupled system being that the smaller component 
of the asymmetric mode demonstrates a split-peak structure in Fig. HJ^b). This feature originates from the anisotropy 
[i.e., the 6>-dependence in the second term in the integrand of Eq. ([T])] of the DDI between the c ores. 

Results of the systematic analysis of the strongly-coupled system are summarized in Fig. 5(a)| and |5(b)[ in the form 



of stability diagrams for all the three types of the modes, in the plane of (a,P), for two different fixed values of h. 
There are two bist ability areas, where the stable flat state and asymmetric solitons coexist [the large yellow triangular 
and trapezoidal regions in the left parts of Figs. |5(a)| and |5(b)[ respectively], or the symmetric and asymmetric 
solitons are simultaneously stable (small orange regions in the right parts if the same figures). The presence of the 
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FIG. 5: (Color online) (a), (b) Stability diagrams in the plane of the total norm (P) and separation a between the cores of the 
strongly-coupled system, for two different values of the overall diameter, h. Symmetric solitons are stable in the bottom right 
areas colored by red, flat states in the bottom left areas, shown in blue, and asymmetric solitons — in the green top areas on the 
left and right sides. The middle-left yellow areas are regions of the bistability area of flat state and asymmetric solitons, while 
the small orange areas at the right center harbor the bistability of the symmetric and asymmetric solitons. (c) The bifurcation 
diagram for solitons in the strongly-coupled system: The soliton's asymmetry, measured by the deviation of the share of the 
total power in one core (Pi/P) from 0.5, versus total norm P. The type of the symmetry-breaking bifurcation is subcritical 
here. The dashed segments, which link the stable asymmetric branches and the stable symmetric one, designate the actually 
missing unstable branches, which, as usual [s^, could not be found by means of the imaginary-time-integration method. 



latter bistability area implies that the SSB transition between the symmetric and asymmetric solitons is subcritical 
in the strongly-coupled system, as explicitly shown in the bifurcation diagram in Fig. [5(H)] [cf. Fig. [3(^. It is thus 
concluded that the repulsive DDI between the cores plays an increasingly more important role with the decrease of the 
separation between them, a, which becomes the dominant parameter of the system at a < (2/3)6. In particular, only 
in this region the symmetric flat state may be stable. On the other hand, the effect of parameter a is inconspicuous 
at a > (2/3)6, like in the weakly-coupled system, see the previous subsection. Accordingly, at a 6, t he are a of the 
bistability between the symmetric and asymmetric solitons (the small orange area in Figs. 5(a)| and |5(b)[ ) shrinks 
to zero, which demonstrates that the SSB changes to the supercritical type, like in the weekly-couple case. Finally, 
antisymmetric solitons and flat states can be also found as stationary solutions of the strongly-coupled system, but 
both these species of the modes turn out to be completely unstable. 



D. Effects of the local nonlinearity 



The inclusion of the local self-attractive {g < 0) or self-repulsive {g > 0) nonlinearity into the coupled GPEs, Eq. 
([1]), shifts the point of the SSB bifurcation (the critical value of the total power. Per), but it does not change the 
supercritical character of the bifurcation in the weakly-coupled system. The bifurcation diagrams for the system of 
this type, in the presence of the local nonlinearity of either sign, are presented in Figs. [6Ka) and (b), cf. Fig. [3l 
In addition. Fig. ^c) displays the dependence of the bifurcation point. Per, on strength g of the local nonlinearity. 
Naturally, the latter dependence is weaker for smaller values of the pipes' diameter, 6, as the DDI is stronger at 
smaller 6, suppressing the effect of the local nonlinearity. It is natural too that the self-attractive local nonlinearity 
{g < 0) makes Per smaller, while the local self-repulsion {g > 0) makes it larger. Similarly, the addition of the contact 
nonlinearity does not change the type of SSB bifurcation in the strongly-coupled system either (not shown here in 
detail). 



III. MOBILITY AND COLLISIONS BETWEEN SOLITONS 



Collisions between solitons represent an important aspect of the dynamics of integrable and nonintegrable models 
[rij . In dual-core systems, one can study collisions between symmetric solitons, and, which is most interesting, 
between asymmetric ones with equal or opposite polarities^ i.e., with larger components belonging to the same or 
different cores [72[. These three cases are schematically defined in Fig. [71 

We here report numerical results for collisions between solitons in the weakly-coupled system, taking the initial 
state at t = as a pair of far-separated kicked solitons, 

^1,2 = /7}^^2\x + xo,P)e^^(^+^°) + /7}^2^(x-xo (7) 
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FIG. 6: (Color online) The bifurcation diagrams in the weakly- coupled system, with (a, 6) = (1,0.4), in the presence of the 
self-attractive (a), ^ = —2, or self-repulsive (b), ^ = +2, local nonlinearity. (c) The total norm at the bifurcation point as a 
function of strength g of the additional contact nonlinearity, at different fixed values of h. 
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FIG. 7: Three types of collisions between solitons in the dual-core system: (a) "Unipolar" asymmetric solitons, with the 
larger components belonging to the same core, (b) Asymmetric solitons with opposite "polarities", (c) The collision between 
symmetric solitons. 



where 2 {x ^ xq^ P) are the two-component soHton solutions with total powers P, which are centered at x = T^o, 
with sufficiently large initial separation 2xo, and 77 is the kick (momentum imparted to the soliton). Solutions /7^^^ 
and [7(2) are either identical ones, or, in the above-mentioned case of opposite polarities, these are two asymmetric 
solitons with swapped components, see Fig. [Tl^b). 

In Figs. [8] and [9l we display typical results of collisions for asymmetric soliton pairs in the case which is close to 
the border between the weakly- and strongly-bound systems, (a, 6, P) = (1,0.4,0.5). In panels [51^ a, b) and[9fa,b), it is 
observed that slowly moving solitons bounce back from each other elastically. When the kick and ensuing velocities 
are larger, the collision becomes inelastic, leading to the merger of the solitons into a single asymmetric breather. 
Panels [8l^c,d) and [8l^e,f) demonstrate, severally, that the merger of the colliding unipolar solitons may switch the 
polarity of the emerging breather, in comparison with the original solitons, or produce a breather which keeps the 
original polarity. The merger of the solitons with opposite polarities gives rise to an asymmetric breather, whose 
polarity is established spontaneously, as seen in Figs. [9l^c,d). Finally, at still larger values of the kick, the moving 
solitons pass through each other, reappearing in an excited form (each one as a moving breather), see Figs. [8fg,li) 
and[9fe,f). The norms of the outgoing vibrating solitons are equal, as they were before the collision. 

The same sequence of outcomes of the collisions — rebound, merger into a breather, and passage in the form of 
moving breathers — is observed, with the increase of the kick (77), at other values of the parameters. Collisions 
between symmetric solitons seem simpler (not shown here in detail): Rebound at small values of 77, and passage, 
without conspicuous excitation of intrinsic vibrations, at larger rj. 



IV. THE PT-SYMMETRIC VERSION OF THE WEAKLY COUPLED SYSTEM 



It was proposed to realize PT-symmetric systems in BEC by linearly coupling two traps (cores, in the present 
terms) with the loss of atoms in one of them, and compensating supply in the other, which may be provided by a 
matter- wave laser [73|- The objective of the analysis was to provide for a direct realization of the PT symmetry in 
quantum media, after it was proposed [7^] and implemented [75| in classical optics, see also review (76| . 
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FIG. 8: (Color online) Examples of collisions between unipolar asymmetric solitons: (a,b) for the slow solitons, with 77 = 0.1; 
(c,d) for the intermediate velocity, 77 = 0.4; (e,f) for a larger intermediate velocity, 77 = 0.8; (g,h) for the fast solitons, 77 = 1.6. 



Accordingly, the T^T-balanced version of hnearly-coupled GPEs ([T]) is 
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where 7 > is the coefficient of the gain and loss in the first and second cores, respectively, = Gdd = 1 is fixed, as 
above, and the local nonlinearity is dropped = 0). As well as in the recently studied model of the T^T-symmetric 
coupler with the cubic nonlinearity , stationary T^T-symmetric and antisymmetric solutions to Eq. (0) can be 
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FIG. 9: (Color online) Examples of collisions between asymmetric solitons with opposite polarities: (a,b) for the slow solitons, 
with ?7 = 0.1; (c,d) for the intermediate velocity, 77 = 0.4; (e,f) for the fast solitons, 77 = 1.6. 



found as 



iIj[^2^^\x ^ t) = ±ze~*^^0(x) exp ^^^-z arcsin7 j , 



(9) 
(10) 

where the upper and lower signs correspond to subscripts 1 and 2, respectively, /i is a real chemical potential, and real 
function is the solution of stationary equation (|2]) for symmetric solitons in the system without the VT terms, 
and with the same value of /i. Obviously, the symmetric and antisymmetric solitons exist for 7 < 1, and they form 
continuous families parameterized by /i, like their counterparts in the conservative system. 

Here we focus on the analysis of the stability of the T^T-symmetric and antisymmetric solitons (|9]), which is a 
nontrivial problem in the present context. Typical examples of stable symmetric solitons are displayed in Fig. [TOl 
along with their counterpart in the system without the VT terms (7 = 0). 

The PT-symmetric and antisymmetric solitons remain stable up to a certain critical value, 7cr, of the gain- loss 
coefficient, and are unstable in the interval of 7cr < 7 < 1, which is a situation typical for solitons in T^T-symmetric 
systems [49|, [50|, [l3 • The unstable solitons suffer a blowup of the pumped component and decay of the damped one, 
which is a typical scenario too (not shown here in detail). The most essential results are presented in Fig. [TTI for 
the weakly-coupled system, in the form of the dependence of 7cr on the total norm for different diameters of the 
parallel-coupled pipes, b (recall b was demonstrated above to be the most essential coefficient for the weakly-coupled 
system). The figure demonstrates that the instability region expands as the nonlinear interactions get stronger, which 
is caused either by the increase of the total norm (P), or decrease of b. 



V. CONCLUSION 



The objective of this work is to extend the study of the symmetry breaking of solitons in dual-core systems with 
the cubic nonlinearity and linear coupling between the cores, that has been previously analyzed in full detail for local 
interactions, to dipolar BEC with the long-range interactions, which act both inside each core and between them. 
Two versions of the system were introduced, weakly- and strongly-coupled ones, depending on the relation between 
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FIG. 10: (Color online) Examples of stable T^T-symmetric solitons, defined as per Eq. (|9]), with (a, 6, P) = (1,0.4,0.3), for 
7 = 0, 0.3, and 0.6. (a) The profile of \ipi,2{x)\, which is common for the three solitons. (b,c) Real and imaginary part of the 
solitons. 
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FIG. 11: (Color online) Critical value 7cr (the stability boundary for the PT-symmetric solitons) versus the total norm, P, 
at different fixed values of the pipe's parameter, 6, in the weakly- coupled system, (a) For the symmetric solitons (|9]); (b) for 
antisymmetric ones (p^ . The continuous and dashed curves are guides for eye. 



the diameter of the pipe-shaped traps and the distance between them. In either case, the linear coupling accounts 
for hopping of atoms between the cores. The symmetry-breaking bifurcation and stability regions for symmetric 
and asymmetric solitons, as well as for non-bifurcating antisymmetric ones, have been identified in both systems. 
In addition to the solitons, the strongly-coupled system supports a stability region for flat states with the unbroken 
symmetry between the cores, due to the competition between the attractive and repulsive intra- and inter-core dipole- 
dipole interactions. Collisions between kicked asymmetric solitons in the weakly-coupled system were systematically 
studied too, showing bouncing back at small velocities, merger into an asymmetric breather in the intermediate range, 
and reappearance of vibrating fast solitons after the collision. Finally, a PT-symmetric generalization of the weakly 
coupled system was introduced as a more "exotic" extension of the system, and a stability boundary for T^T-symmetric 
and antisymmetric solitons was found. 

A challenging generalization of the analysis may be its application to 2D dual-core systems, where, in particular, 
the symmetry breaking may occur not only for fundamental solitons, but also for solitary vortices. 
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